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Abstract 

The quantum ratchet effect in fully chaotic systems is approached by studying, for the first time, 
statistical properties of the ratchet current over well-defined sets of initial states. Natural initial 
states in a semiclassical regime are those that are phase-space uniform with the maximal possible 
resolution of one Planck cell. General arguments in this regime, for quantum-resonance values of 
a scaled Planck constant h, predict that the distribution of the current over all such states is a 
zero-mean Gaussian with variance ~ Dh? /(2-7T 2 ), where D is the chaotic-diffusion coefficient. This 
prediction is well supported by extensive numerical evidence. The average strength of the effect, 
measured by the variance above, is significantly larger than that for the usual momentum states 
and other states. Such strong effects should be experimentally observable. 
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I. INTRODUCTION 



Understanding quantum transport in generic Hamiltonian systems, which are classi- 
cally nonintegrable and exhibit chaos, is a problem of both fundamental and practical 



importance. The study of simp 
quantum-transport phenomena 



experiments |2 , 



3 
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e model systems have led to the discovery of a variety of 
lMl8l|. several of which have been observed in atom-optics 



- |l8| and allow to control the quantum motion of cold atoms or 



Bose-Einstein condensates in different ways. Recently, classical and quantum Hamiltonian 

- Il5| and 



ratchets" have started to attract a considerable interest, both theoretically 10 
experimentally 16l-ll8|. Ratchets are usually conceived as spatially periodic systems with 
noise and dissipation in which a directed current of particles can emerge from an unbiased 
(zero- mean) external force due to some spatial/temporal asymmetry 19] . In classical 
Hamiltonian ratchets [10], dissipation is absent and noise is replaced by deterministic 
chaos. A basic general result for Hamiltonian dynamics under an unbiased force is that 
the average current of an initially uniform ensemble of particles in phase space is zero 10 ]. 
As a consequence, a completely chaotic system carries essentially no ratchet current. On 
;he other hand, the corresponding quantized system can feature significant ratchet effects 
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18| . An important problem is to understand the nature of these full-chaos quantum 



effects in a semiclassical regime, in particular how precisely they vanish, as expected, in 
the classical limit. All the studies of quantum chaotic ratchets until now have mainly 
focused on the impact of several kinds of asymmetries on the quantum directed current 



from a fixed initial state, 
the initial state [lO 
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t is, however, well established that the current is sensitive to 
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171 . Il8| and this sensitivity is expected to be especially high 



in a semiclassical full-chaos regime, reflecting the exponential sensitivity of chaotic motion 
to initial conditions. Thus, to get a comprehensive understanding of the quantum ratchet 
effect, it is necessary to adopt a more global approach, not limited to a single initial state. 



In this paper, the semiclassical full-chaos regime of quantum ratchets is approached by 
studying, for the first time, statistical properties of the current over sets of initial states 
with well-defined natural characteristics. The systems considered are generalizations of the 
paradigmatic kicked Harper models (KHMs) 



14 



15], with Hamiltonian 



H = L cos(p) + KV(x) 6 ^ ~ *)> 



(1) 



where L and K are parameters, x and p are scaled position and momentum operators, V(x) 
is a general 27r-periodic potential, t' is time, and t is the integer time labeling the kicks. 
Generalized KHMs such as ([T]) describe several realistic systems , |l5j], in particular 



they are exactly related 6|, |7| to kicked harmonic oscillators, which are experimentally 
realizable by atom-optics methods 3J], and to kicked charges in a magnetic field [7]. 
Recently 14|, the systems ([1]) were shown to exhibit generically a significant and robust 
quantum momentum current (ratchet acceleration) under full-chaos conditions. The initial 
state was chosen, as in other works, as a zero-momentum state. In our statistical approach, 
we identify natural initial states for the semiclassical regime as those that are analogous as 
much as possible to a phase-space uniform ensemble, for which classical ratchet effects are 
totally absent. These are states which are uniform in phase space with the maximal possible 
resolution of one Planck cell. Such uniformity is not featured by a momentum state which 
is uniform in position but is infinitely localized in momentum. 



Assuming quantum-resonance values of a scaled Planck constant H — [x,p]/i in the semi- 
classical regime, we derive an estimate for the distribution of the quantum momentum 
current / over maximally uniform initial states: This distribution is a Gaussian with mean 
(J) = and variance (A/) 2 = (I 2 ) ~ Dfr 2 / (2tt 2 ), where D is the chaotic-diffusion coefficient. 
A good agreement is found between this estimate and extensive numerical results using an 
exact formula for / which we also derive. Examples are shown in Fig. 1 and will be discussed 
further in Sec. IV. The average strength of the effect, measured by the variance above, is 
found to be significantly larger than that for the usual momentum states and other states 
exhibiting also zero-mean Gaussian current distributions (see the insets of Fig. 1 and Sec. 
IV). Our results should be experimentally observable to some extent using states approxi- 
mating the maximally uniform states. 



The paper is organized as follows. In Sec. II, we define maximally uniform states in 
phase space. The main result, i.e., an estimate of the momentum-current distribution 
over these states in a semiclassical full-chaos regime, is derived in Sec. III. Numerical 
evidence for this result is provided in Sec. IV, where we also study the momentum-current 
distributions of states which approximate the maximally uniform states; momentum states 
are the crudest approximating states. Conclusions are presented in Sec. V, where we briefly 
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I/AI 



FIG. 1: (Color online) Distributions of the normalized quantum momentum current I/AI over 
maximally uniform states for h = 2-7r/121 in two extreme cases of fully chaotic systems ([I]) with 
K = 15: The symmetric case 11 S" with V{x) = cos(x) and L = K (red squares, AI = 0.086) and 
the strongly asymmetric case "A" with V{x) = cos(x) + sin(2x) and L = K/2 (blue diamonds, 



l4 | for a zero-momentum initial state. The 



AI = 0.214); the latter case was studied in Ref. 
origin of ratchet currents in case S is explained in Sec. IV. The insets show the distribution of 
I/AI over momentum states in case A with K = 20 (left inset, AI = 0.043) and over low-order 
approximations of the maximally uniform states in case S with K = 15 (right inset, AI = 0.026), 
see Sec. IV for more details. The solid line in all plots is a zero-mean Gaussian with variance 1. 

mention possible experimental realizations of the strong quantum-ratchet effects predicted. 
Detailed derivations of some exact results are given in the Appendix. 



II. MAXIMALLY UNIFORM STATES IN PHASE SPACE 

Maximally uniform states in phase space are denned on the basis of the translation opera- 
tors T x (a) = exp(ipa/h) and T p (b) = exp(—ixb/h) shifting x andp by a and b, respectively. A 
state is uniform on the phase-space lattice with unit cell formed by a and b if it is invariant 
under application of T x (a) and T p (b) up to constant phase factors: T x (a) = expect) 
and T p (b) = exp(— i(5) \ip). This means that \ip) is a simultaneous eigenstate of T x (a) and 
T p (b), so that these operators must commute. Using T x (a)T p {b) = exp(—iab/h)T p (b)T x (a), 
we see that [T x (a),T p (b)] = only if ab is a multiple of h = 2nh. Maximally uniform states 
|^) correspond to the smallest unit cell, i.e., the Planck cell with area ab = h. In this case, 
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one can easily check that the position representation of the eigenstates is explicitly given 



by (20| 

oo 

(x |Vw> = -7= Yl e 2 ™ /b 5(x - Wl - na), (2) 

v n=— oo 

where w = (u>i,io 2 ) range in the Planck cell, < W\ < a, < w 2 < b, and specify the 
phases a and (3 above: a = W2a/h and = w\b/h. The states (T5]) for all w form a complete 



and orthonormal set 



20(| . Simple choices of (a, b) can be made by observing that the 



one-period evolution operator for (JTJ, [7 = exp[— L cos (p) /ft] exp[—KV(x)/h], is 27r-periodic 
in both (x,p). Thus, the torus T 2 : < x,p < 2n is a reduced phase space for the system. 
For simplicity, we shall assume from now on that there are precisely an integer number N 
of Planck cells within T 2 , choosing a = 2tt/N and b = 2tt. Then, h = ab/(2n) = 27r/N, 
so that the semi classical regime S<1 corresponds to N 3> 1. As it is well known {4, 8], 
the values 2tt/N of h are those for which a classical-quantum correspondence can be most 
easily established for systems describable on a phase-space torus. These values correspond 
to the main quantum resonances in the semiclassical regime. 



III. SEMICLASSICAL ESTIMATE OF 
THE MOMENTUM-CURRENT DISTRIBUTION 

The momentum-current operator I can be formally defined in the Heisenberg picture as 
I = lim^oo ptj 1 /t, for integer time t. Then, the momentum current /(w) for initial state 
((2j) is the expectation value of I in (j2j). More precisely, we show in the Appendix that in the 
basis of states (j2J) / is essentially represented by (-?/>w|-f l^/v) = I(w)S(wi — w' 1 )5(w2 — w'2), 
where /(w) is given by the explicit exact formula: 

/(w) = -»£>(0iw)|'?§£=>. (3) 

i=i 1 

Here 4>j(0; w), j = 1,...,A, are coefficients appearing in expressions connecting states 
02]) with the iV quasienergy (Floquet) eigenstates |^j, w ) of the evolution operator U for 
= 2n/N and Ej(w) are the corresponding quasienergy levels. 
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We now derive from general arguments the following estimate for the distribution T(I) 
of /(w) over w in a semiclassical full-chaos regime: 



r(J) ~ ; exp 

v^AJ 



. rx2 2D Dh 2 
( A/ ) ~ ThT = "o3T> ( 4 ) 



2(A/) 2 J ' v 7 iV 2 2vr 

where D is the chaotic-diffusion coefficient. To derive @, we first identify natural classi- 
cal analogs of the states fl2]). To this end, let us calculate the momentum representation 
{p l^w) = exp(— ipx/h)(x IV'w) dx of \if) w ). One has, up to an irrelevant constant factor, 

oo 

(P l^w) = Yl e- 2mnWl/a S(p -w 2 - nb). (5) 

n=— oo 

It is clear from the delta combs (T5]) and (j3J) that |-0 W ) is associated with the phase-space lattice 
(x,p) = w + z(n), where z(n) = (rixa, n 2 b) = (27rni/N,2im 2 ) for all integers n = (ni,^). 
In fact, one can easily show that the Husimi distribution of fl2]) is peaked on every point 
of the lattice w + z(n). This lattice, viewed as an initial phase-space ensemble, is thus a 
classical analogue of |t/> w )- Next, consider the classical one-period map M for the systems 
(JTJ) : pt+i = Pt+Kf(x t ), x t +i = x t — Lsm(p t+ i), where f(x) = —dV/dx. For initial conditions 
z o — ( x o,Po), the classical momentum current in t iterations is I c ,t( z o) — Ap t (z )/t, where 
Ap t (z Q ) = p t —po. Since the map M is clearly 27r-periodic in both x and p, one can restrict 
z to the phase-space torus T 2 : < x,p < 2tt. Accordingly, the initial ensemble w + z(n) 
with z(n) = (27rni / 'N , li^n^) can be restricted to a finite lattice of N points in T 2 with 
rii = 0, . . . , N — 1 and n 2 = 0. The classical analog of the quantum current /(w) is the 
average I c ,t( w ) °f ^c,t(zo) = Ap t (z )/t, with z = w + z(n), over this finite lattice: 

^(w) = ^^A Pt [w + z(n)], (6) 

n 

for some time t to be specified below. Now, under strong-chaos conditions (large K and 
L) and for sufficiently large t, each of the iV quantities Ap t [w + z(n)] in (jBJ) is expected to 
behave diffusively, i.e., to be distributed over w approximately as a Gaussian with mean 
(Ap t ) = and variance ((Ap t ) 2 ) ~ 2Dt. Since these N quantities are associated with iV 
different initial points w + z(n) in the chaotic region, they should behave as independent 
(uncorrelated) random variables. It then follows from the central limit theorem that for 
large enough iV the average current (jSJ) is distributed over w as a Gaussian with (l c ,t) = 
and (I 2 t ) w N ((Apt) 2 ) /{Nt) 2 w 2D/(Nt). This shows how (l 2 tt ) decays to zero as t -t oo, 
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when chaotic orbits explore ergodically all the phase space. Since there are precisely iV 
Planck cells in T 2 , a typical such orbit will explore phase space, after a time t ~ N, up to 
the maximal quantum resolution of one Planck cell. Then, the distribution of /(w) over w 
is expected to be approximately the same as that of the classical currents for t = N, 
i.e., a zero-mean Gaussian with variance (A/) 2 = (I 2 ) ~ 2D/N 2 ; this is Eq. (jlj). 



IV. NUMERICAL EVIDENCE AND APPROXIMATING STATES 

In this section, we provide numerical evidence for the semiclassical estimate (jl]) using 
the exact formula ([3]) and study the momentum-current distributions for states which 
approximate the maximally uniform states; at the lowest order of approximation, the 
approximating states are just momentum states. First, the distribution T(I) was calculated 
using ([3]) for several potentials V(x) and many large values of K and L in A good 
agreement was generally found between T(I) and a zero- mean Gaussian distribution for 
sufficiently small h. As representative examples, Fig. 1 shows distributions of I/AI in two 
extreme cases U S" and U A" described in the caption. For the assumed quantum-resonance 
values 2tc/N of h, the origin of ratchet currents in the symmetric case S is the same 



as that already established in recent theoretical [13J and experimental 13, ll8] works on 
quantum-resonance ratchets: This is a relative asymmetry caused by the non-coincidence 
of the symmetry centers of a symmetric potential with those of a symmetric initial state. 
The potential V(x) = cos(a;) in case 5* has symmetry centers at x = 0, 7r while the state ([2]) 
has them at x — Wi, w± + a/2. Thus, for generic values of Wi, /(w) ^ 0. 

Natural approximations of the maximally uniform states, denoted in what follows by 
ip^l for integer B, correspond to truncations of the infinite sum in ([5]): 

B 

(P } > = E e- 2m/a 5(p -W2- nb). (7) 

n=-B 

The states (j7j) are superpositions of the 2B + 1 momentum states \p = W2 + nb), |n| < B, 
and should be experimentally realizable (see next section). In particular, ip™^ are j us t 
momentum states with p = w<i- We denote by 7b(w) the momentum current for initial 
state ([7j) and by (AIb) 2 the corresponding variance over w. The left inset of Fig. 1 shows 
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the distribution of Jb(w)/A/b over w for B = (momentum states) in case A while the 
right inset shows it for B = 2 in case S; the currents for momentum states vanish in case 
S. The fact that these distributions are again approximately zero-mean Gaussians could be 
expected from the simple relation ( JT5"j) between 7b(w) and /(w) derived in the Appendix. 
Actually, we show in the Appendix that this relation leads to the exact result 

(a/b) 2 = iji) < (Air = (n ■ (8) 

The quantities AI and AIb were extensively studied as functions of several pa- 
rameters. Consider the naturally normalized variance R = N 2 (AI) 2 /(2D q \), where 
D q i = K 2 [f (x)] 2 dx / 2 is the "quasilinear" value of the diffusion coefficient D, obtained 
from the KHM map M above by neglecting all the force-force correlations C t = (f{xo)f{xt)), 
t 7^ 0; for sufficiently strong chaos, D is very close to D q \. The semiclassical estimate for 
the variance in Eq. (jlj) would imply that R ~ D/D q \. Indeed, Fig. 2 shows a reasonably 
good agreement between R and D/D q \ versus K in both cases 5* and A. Discrepancies 
arise mainly around peaks of D/D q \, especially the peak near K 6.5 in case S, due to a 
small accelerator-mode island. Thus, for general large K with R rj 1, AI increases almost 
linearly with K like ^jD q \. 

Fig. 3 shows loglog plots of AI versus iV in cases S and A. The results agree very 
well with the N~ l behavior predicted by (Jl|). Fig. 4 shows plots of AIb/ AI versus B in 
the two cases. We see that AIb is always smaller than AI, in accordance with the exact 
inequality OH]), and approaches monotonically AI as the order B of approximation increases. 
For small B, AIb is significantly smaller than AI and attains its minimal value at B = 0, 
corresponding to momentum states (A/o = in case S). Within the limited domain of iV 
we could study numerically, AIb appears to decay with N like iV M , where /i ranges between 
—0.9 to —1.1 with an error not smaller than ±0.03. 

V. CONCLUSIONS 

In conclusion, we have presented a first study of the semiclassical full-chaos regime of the 
quantum ratchet effect using a novel statistical approach which is most required in view of 
the sensitivity of the effect to the initial state. For maximally uniform states l^w) m phase 
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FIG. 2: Filled circles: The quantity i? = Af 2 (A/) 2 /(2-D q i) versus K for iV = 2vr/?i = 121 in cases 
S [(a), D q i = K 2 /4] and A [(b), D q i = 5/f 2 /4] denned in the caption of Fig. 1. Crosses joined by 
a line: D/D q \ versus K. 
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FIG. 3: (Color online) Loglog plots of AI versus N = lixjh (49 < N < 169) for K = 15 in case 
S (red squares, with linear fit having slope fj, = —0.98 ± 0.01) and in case A (blue diamonds, 
H = -1.02 ±0.01). 



space, used here for the first time as natural initial states, the momentum-current distribu- 
tion (jlj) exhibits clear fingerprints of classical chaotic diffusion, a genuine quantum-chaos 
phenomenon. The simple formula for the variance in Eq. (jlj) involves just D and h 2 and is 
thus interestingly similar to the well-known one for the localization length £ in the kicked 
rotor 2l[, £ ~ D/h 2 (in our notation). This variance was shown to be significantly larger 



than that for momentum states (B = in Fig. 4), which were standardly used in previous 
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FIG. 4: (Color online) AI^/AI versus B for iV = 121 and K = 15 in case S (red squares) and in 
case A (blue diamonds). 

works within the ordinary (single-state) approach. The states then appear to give the 
strongest quantum ratchet effect known until now. 

One can approximate |^ w ) to order B by the states ([7]) which are superpositions of 2B + 1 
momentum states and whose variance increases with B (see Fig. 4). Superpositions of two 
momentum states were recently used in experimental realizations of quantum-resonance 



ratchets 
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181 ] . The states ([7]) can be experimentally prepared for at least B < 10 22]. 
Also, the systems (CE]) are exactly related p, 0] to kicked harmonic oscillators which are 
experimentally realizable {3)]. Thus, the current distributions and the strong quantum 
ratchet effects predicted in this work should be observable to some extent in the laboratory. 
More detailed aspects of our statistical approach and its extension to other systems and 
parameter regimes will be considered in future studies. 
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APPENDIX 

We derive here formula fl3]) and the inequality (|8]). We start with a summary of results 
from Ref. js]. Since the evolution operator U = exp[— L cos(p)/H] exp[— KV(x)/h] is 2tt- 
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periodic in (x,p), it commutes with both T x (2ir) = T^(a) (a = 2n/N) and T p (b) (b = 2%). 
Therefore, one can find simultaneous eigenstates of U, T^(a), and T p (b). The general 
eigenstates of T^(a) and T p (b) are given in terms of (J2]) by 

N-l 

l*J.w) = ^ <M m ; W ) IVV,™2+ma) • (9) 
m=0 

Here 0j(m;w), j = 1,...,JV, form A 7 independent vectors of coefficients, 
Vj(w) = {0j(m; w)} m ~Q, which are determined from the eigenvalue equation 
U l^j.w) = exp[— iEj(w)] |^j, w ), where Ej(w) are the quasienergies. It is clear from Eq. (J9]) 
that Vj(w) are the representation of l^w) i n the iV-basis \4 , Wlt w 2 +ma) , m = 0, ...,N — 1. 
ii this basis, U is represented by an N x N unitary matrix M(w) with known elements 
. Thus, M(w)Vj(w) = exp[—iEj(w)\Vj(w). This completes the summary of relevant 
results from Ref. [81. 



Let us now calculate the matrix element ^-?/> w |/|?/vy, where I = lim^oo U 1 pU 1 /t is the 
momentum-current operator. First, we have 

(^icHpC/VtlVv) = (VwIpIVv,*) A, (io) 

where |Vw) = ^* IV'w}- Using the completeness of the eigenvectors Vj(w) of M(w), i.e., 
J2f=i <^( m ; w)(f) j (m'; w) = 5 m . m ,, we can invert Eq. © to get \ip w ) = J2f=i w ) l*i,w)- 
One then has 

N 

To determine the asymptotic behavior of Eq. ffTT]) for large £, we use the equation 
tf t {x : p) |^/j, w ) = exp[— 2iEj(w)] |^j, w ), the expansion (Q, and the fact that p = —ihd/dx is 
represented by —ihd/dwi in the basis (jSJ) [20], due to the delta comb in x. Using also the 
orthonormality of ([2j), (VvlVv) = — w[)5(w2 — w 2 ) and of Vj(w), we find that 
the dominant terms in ( JJ t '^>ji vv \p\ tj t ^!j^i \ for large t give the asymptotic behavior 



f>%, iW |p| £/%, w ') ~ -th dE ^ Sjj.Siwx - w[)5(w 2 - w' 2 ), t > 1, (12) 

where we assume for simplicity that < W2,w' 2 < a. After inserting ( TT2]) in (TTTj) and 
dividing by t, we get from ( flOl) in the limit £ — > oo: /^wl ^|Vv^ = I(w)5(wi — w' 1 )5(w 2 — w 2 ), 
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where I(w) is given by formula ([3]). 



To derive the inequality (|S]), we first notice that the states §5§ and (J7J) can be easily 
related as follows: 

\^ ] ) = I dw' l g B {w' 1 - wi) IV^.tua) , ( 13 ) 
Jo 

where 

1 B 

ffs(Wl) = 77=== ; = eX P /«) • ( 14 ) 



n=-B 



The factor before the sum in (fl4"l) assures the normalization 

dw^Kwi) = 1. (15) 



o 



Eq. (fl5l) implies, because of (VvlVv) = — u^oXu^ — w 2 ) an d (V'w^lV'v 

I(w)5(wi— w' 1 )5(w2— w' 2 ) (see above), that the states f|T3|) satisfy the orthonormality relation 
(V'StalV'Si^/ = 5 ( W2 ~ ^2) and (V©**!*!^^) = Ib(w)5(w 2 - w 2 ), where 

■Jb(w) = / dw' 1 g 2 B {w[ - twi)/(twi,«;2)- (16) 
Then, (fT6j) is clearly the momentum current for initial state ([TBI . 



< / dw' 1 \F(w' l )\ 2 / d^|GK)f 
Jo Jo 



Now, using f|T5|) . (fT6l) . and the Cauchy-Schwarz inequality 
dw'^iw'^Giw'^ 

with the identifications F(w[) = gB(w[ — Wi)I(w[,W2) and G(w' 1 ) = gB(w[ — Wi), we get: 

7l(w)< I" dw' 1 g 2 B (w' 1 -w 1 )I 2 (w' 1 ,w 2 ). (17) 

Using (fT7|) and, again, (JTB1) in the definition = j Q a dw\ dw2l B (w)/h, with (J 2 ) 

similarly defined, we finally obtain that (J|) < (J 2 ). This is inequality (jSJ), where 
{AI B ) 2 = (1%) since (Jb) = 0, as easily implied by (TIB]) and (I) = 0. 
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